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Abstract. 

We investigate the structure of topological defects in the ground states of 
spinor Bose-Einstein condensates with spin F = 1 or F = 2. The type and 
number of defects are determined by calculating the first and second homotopy 
groups of the order-parameter space. The order-parameter space is identified with 
a set of degenerate ground state spinors. Because the structure of the ground 
state depends on whether or not there is an external magnetic field applied to the 
system, defects are sensitive to the magnetic field. We study both cases and find 
that the defects in zero and non-zero field are different. 



1. Introduction 

Bose-Einstein condensates (BECs) of alkali atoms have an internal degree of freedom 
due to the hyperfine spin of these atoms. If a BEC is realized in a magnetic trap this 
degree of freedom is frozen and in a mean-field limit the condensate is described by a 
scalar order- parameter. However, if an optical trap [Q is used to confine condensate 
atoms, this degree of freedom is liberated and has to be taken into account Ej- 
Condensates with this property are called spinor or vector condensates. In the mean- 
field theory the ground state of a spinor condensate is described by an order-parameter 
\P(r) = y/n(r)£(r), where n(r) is the density of the condensate and £(r) is a normalized 
spinor, £'(r)£(r) = 1. In this paper the density n is assumed to be constant. Because 
of the vectorial nature of the order-parameter, the behaviour of spinor condensates is 
in many ways different from that of scalar condensates. One manifestation of this can 
be seen in the difference of defects in scalar and spinor condensates. In the former 
vortices with integer winding numbers can exist. The latter allows for more complex 
defects, which are the topic of this paper. Our study is based on the ground states 
calculated using the mean-field theory and single condensate approximation OH El El- 
Mean-field theory is widely used in the study of Bose-Einstein condensates and it is 
usually assumed to give a good description of the physical system. However, some 
results suggest that the actual ground states of spinor condensates may be different 
from those obtained using mean-field theory |18l llfjj . Thus, the results of this paper 
are valid only as long as mean-field theory can be applied. 
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2. Characterization of the used techniques 

2.1. The order-parameter space 

In condensed matter systems the concept of an order-parameter is very important 
01 0. Order-parameter /(r) is a continuous mapping from some region of the physical 
space into the order-parameter space M, which consists of all possible values of 
the order-parameter. It is usually possible to associate the order-parameter space 
with a group G that acts on that space. If this action is transitive (i.e for every 
x,x' € M there exists some g <E G for which x' — g ■ x), we can arbitrarily choose 
some element x re f € M which we call the reference order-parameter. Every element 
x € M can then be obtained from x re f by acting on it by a suitable element of 
the group G. Those elements of G which leave £ re f fixed constitute a subgroup 
H called the isotropy group. Explicitly H = {g £ G\g ■ x re f = £ re f}- Under 
some rather general requirements for G and M the order-parameter space M can 
be identified with the quotient space G/H. When considering the defects of spin-F 
Bose-Einstein condensate the order-parameter is the normalized (2F + l)-component 
spinor £(r) € C 2F+1 . Because we study what kind of defects can exist in the 
ground state of the system, the order-parameter space is the set of spinors that 
minimize the energy. In the absence of an external magnetic field we can often choose 
G = U(l) X 50(3). However, this choice is not always the correct one, as there 
may be order-parameter spaces in which U(l) x 5*0(3) does not act transitively; see 
below. 17(1) x 50(3) acts on C 2F+1 via equation ((c,i?),£) i-> c/j( F )(i?)£, where 
(c, R) G U(l) x 50(3) and is the 2F + 1 dimensional irreducible representation 
of 50(3). This representation is given by the map R(a,/3,j) i— > D^ s \a, /3,~f), where 
i?(a,/3,7) = R z (a)R y ((3)R z (j) S 50(3) is given as a product of rotations about y- 
and z-axes and (a, /?, 7) = exp(— iaF z ) exp(— i(3F y ) exp(— i"/F z ). Here F y and F z 
are y- and z- components of the spin matrices corresponding to spin F. Representation 
matrices for F = 1 and F = 2 are given in the appendix. 

In the presence of an external magnetic field the energy of the ground state is 
invariant under gauge transformations and rotations about the axis of the magnetic 
field, so we choose G = 17(1) x 50(2). 

2.2. Homotopy groups 

Homotopy groups of the order-parameter space describe physical defects @j. The nth 
homotopy group 7r n (M) of the space M consists of the equivalence classes of continuous 
maps from n-dimensional sphere 5" to the space M. Two maps are equivalent if they 
are homotopic to one another. In physics, the first and second homotopy groups 
are of special importance. The first homotopy group ni(M) describes singular line 
defects and domain walls, which are non-singular defects. The second homotopy group 
7ra(M) describes singular point defects and non-singular line defects. Thus, identifying 
M with G/H, we can learn much from the possible defects in a physical system if we 
know Ttx(G/H) and ^(G/H). These can be calculated with the help of the following 
theorem. 

Theorem 1 Let G be a Lie group with ttq(G) = tti(G) = ^(G) = 0. Here denotes 
a one element group. Let H C G be a closed subgroup, and Hq C H the connected 
component of the identity. There are isomorphisms 

iri(G/H)SiH/H (1) 
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Table 1. The reference spinors and their general forms for F = 1 spinor 
condensate when the external magnetic field is zero. 





s ref 
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af 


(0,1,0) 


e ie (-e^ la 4= sin /3, cos /3, e ia 4= sin /3) 



and 

7r 2 (G/H)^7r 1 (H ). (2) 

For a proof, see or [15] . 

We cannot use this theorem if G is U(l) x 50(3), because iti(U(1) x 50(3)) = 
ZxZ 2 . This problem can be solved by using R x SU(2) instead of U(l) x 50(3), since 
this group fulfils the requirements of the theorem. The former is a covering group of 
the latter, the covering projection P : R x SU(2) — > t/(l) x 50(3) being given by 
(x,U(a, (3. 7)) I— > (e lx , R(a, (3, 7)) , where i£l and 

Every matrix in 5(7(2) can be written in this form. Sufficient intervals for a, (3 and 7 
are [0, 2tt], [0, 7t] and [0,47r], respectively. 

3. Spin 1 

The ground state structure for F = 1 condensate was calculated by Ho j2] and by 
Ohmi and Machida If the external magnetic field is non-zero the spin-dependent 
part in the energy is £(£) = c(F)| — p(F z )^ + q(F z }^ [JJ]. Here the kinetic energy 
term is neglected in the Thomas-Fermi approximation and (F)^ = £^F£. The constant 
c depends on scattering lengths and the density n, whereas p describes linear and q 
quadratic Zeeman interaction with the external magnetic field. The external field is 
assumed to be directed along the z-axis. 

3.1. Zero external field 

The equation for energy is obtained by setting p = q = 0. Depending on the sign 
of c energy is minimized either by (F)| = 1 or (F)^ = 0. The former is called 
the ferromagnetic (/) and the latter the antiferromagnetic (af) phase. The order- 
parameter spaces corresponding to these phases are = {£ e C 3 |(F)| = 1, £t£ = 1} 
and Af a/ = £ C 3 |(F) 5 = 0, £t£ = 1}. it i s easy to see that U{1) x 50(3) 
acts transitively on these sets. The reference order-parameters and general order- 
parameters obtained from these by a rotation and gauge transformation are shown in 
table [U 

3.1.1. Ferromagnetic phase From table 1 we see that we do not need the angle 
6 because 7 can produce all possible gauge transformations. This means that 
instead of R x SU(2) we can use only SU(2). To find the elements of the isotropy 
group we set ^(a, f3, 7, 0) = ^ c- This gives = {I, —I}. Using theorem 
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Owe get TTi(G/H :f ) = {±1} and 7r 2 (G/-ff^) = 0. The order-parameter space is 
5£/(2)/{±l} = 50(3). This order-parameter space has also been encountered in 
3 He-.4 [12]. 

Physically, these results mean that we can have one non-trivial singular vortex 
but we cannot have any non-trivial monopoles, i.e. singular point defects. The non- 
trivial vortex is a defect in which the overall phase of the spinor changes by 2ir as the 
defect line is encircled. 

3.1.2. Antiferromagnetic phase Now the isotropy group is H a f = {(n2ir, a((p)), ((n+ 
ga((p)) | (f £ [0,47r],n £ Z}, where we have defined a(<p) = U(tp, 0,0) and 
g = U(0,ir,0). The connected component of the identity is Hq = {(0,a(<p)) \ tp £ 
[0,4tt]}. We get m(G/H a f) S {(n27r, l)ff Q a/ , {{n+\)2-K,g)H% f | n € Z}. This group is 
isomorphic to Z, the isomorphism being given by the map ((2n+ j)7r, g^H^ i— » 2n+j, 
where j — or 1 and g° = I. Thus tti (G/H a f) = Z. Previously the order- 
parameter space and first homotopy group were concluded to be U(l) x 5 2 and Z [2] or 
[/7(1) x 5 2 ]/Z2 and Z x Z2 ;6J. Both of these order-parameter spaces are incorrect but 
the first homotopy group in [5] is correct. However, that of reference P| is not correct, 
since there is no (group) isomorphism between Z and Z x Z2. The isotropy group (in 
U(l) x 50(3)) is isomorphic to 0(2), but it can not be expressed as a direct product 
of a subgroup of U(l) and 50(3). Thus the order-parameter space can be written 
only as G/H = \U(l) x 50(3)]/0(2)g+s, where G + S means that the isotropy group 
consists of gauge transformations performed simultaneously with spin rotations. 

Because is homeomorphic to U(l) and homeomorphic spaces have the same 
homotopy groups, we get 7r 2 (G/7J a ^) = 7r 2 (J7(l)) = Z. If we move around a closed 
path in the condensate we note that when we return to the starting point the angle 
9 has changed by some amount. If we define the change in this angle divided by 2n 
to be the winding number, we see from the elements of H a f / that the winding 
number can be either an integer (n) or a half- integer {\ + n) Paths in the 

order-parameter space can be represented pictorially as follows. From table 1 we see 
that we have three parameters in the general expression for the ferromagnetic state. 
From these a and (3 can be restricted to the intervals [0, 2ir] and [0,7r], respectively, 
and 9 £ [0,2tt]. However, because £ a ^(a ± it, ir — (3, 7, 9) = — t; a * (a, /3, 7, 9), we can 
actually restrict 9 to the interval [0, ir}. These parameters can be represented using 
cylindrical coordinates (a,r,z), where now r = f3, z = 9, see figure 

In summary, possible line defects are those in which the overall phase changes by 
2mi as the defect line is encircled and those in which a phase change of n + 2-Kn is 
accompanied by a 180° spinor rotation. Also point defects, labelled by integers, are 
possible. 

3.2. Non-zero external field 

Straightforward minimization of energy gives four ground states which are degenerate 
with respect to one or two phase variables, see table The identification of the 
order-parameter space G/H is easier than in the absence of the magnetic field. 

In /i,/2 and af states G — U(l) and H = 1, so G/H = U(l) for which 
7Ti([/(l)) = Z and 7T2(t/(l)) = 0. Thus we can have singular vortices with arbitrary 
integer winding numbers but we do not have singular point defects. This resembles 
the situation in a scalar condensate, where we have similar defects. 
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Figure 1. Pictorial representation of the parameters of the antiferromagnetic 
spinor. Angles a and 9 are shown in the picture and angle is the length of the 
vector. On the boundary points (a, /3, 0) and (a ± n, tt — j3, it) correspond to the 
same value of the order-parameter. This is also true for all points for which = it, 
9 is fixed and a £ [0, 2ir). A(B) denotes the starting (ending) point of the curve. 
We assume £(A) = £(B), so the curves define closed curves in the order-parameter 
space. Defining the direction of increasing 9 to be positive, we see in (a) a path 
with winding number | and in (b) a path with winding number —1. 



Table 2. The ground states for F = 1 spinor condensate when the external 
magnetic field is non-zero. States are degenerate with respect to angles 9 and 4>. 
We have assumed that c ^ 0. fs evolves to /i (/2) state as p reaches 2c (—2c). 
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In / 3 state G/H = U{1) x U(l) for which ir x {U{l) x 17(1)) Z x Z and 
7T2([/(1) x = 0. Now we can have independent vortices in the m = 1 and 

m = —1 components of the spinor. 

4. Spin 2 

The ground states for F = 2 spinor condensate were calculated by Ciobanu et al 
and Ueda and Koashi [§]. In the Thomas-Fermi approximation the spin-dependent 
energy is given by £(£) = c(F)| + d|©{| 2 — p(Fz)z, where c and d are constants 
depending on scattering lengths and p describes the linear Zeeman effect. The possible 
ground states are characterized by two parameters, namely ||(F)^|| = ((F)!) 1 / 2 and 
|e c | = |26£-2-2^i+£ 2 |- 

4-1- Zero external field 

The energy in zero magnetic field is obtained by settings = 0. Because (F) 2 and |6^| 
are invariant under the action of U(l) x 5*0(3), and for each £ there exists a rotation R 
for which (F) 2 - = {F Z ) 2 D ^^, we can write the energy in the form £(£) = c{F z )^+d\Q^\ 2 . 
This equation has been solved in 0|B]. It turns out that there are three possible phases 
in the system, ferromagnetic (F,F'), cyclic (C) and polar (P). 

To express the order-parameter spaces we define M(i,j) = {£ € C 5 | ||(F)^|| = 
i, \9^\ — j, £T£ = 1}. Then in F phase the order-parameter space is M(2, 0), in F' phase 



Topological defects in spinor condensates 



6 



Table 3. Ground state spinors and their energies of F = 2 condensate when 
the external magnetic field is absent. General forms of the ground states can be 
obtained from these by a rotation and a gauge transformation. Notice that in P 
state we have two free parameters. Those are needed because we cannot obtain 
every possible spinor representing the polar state from a fixed reference spinor by 
a rotation and a gauge transformation. 



F (1,0,0,0,0) 4c 

F' (0,1,0,0,0) c 

C i(l, 0,^,0,-1) 

P -j= I sin (j> sin?/), sin</> cos tp, vl cos <f>, — sin </>cos ip, sin (/jsini/) ) d 



M(1,0), in C phase M(0, 0) and in P phase M(0, 1). Representative spinors from 
these sets and their energies are shown in tabled It turns out that U(l) x 50(3) acts 
transitively on the parameter spaces of the ferromagnetic and cyclic phases. However, 
this is not true for the order-parameter space of the polar phase. For example, if we 
first choose (j) — and then <j> = = in the spinor representing a polar state, 
we get two spinors which can not be converted to each other by a rotation and gauge 
transformation. This means U(l) x SO(3) is not a group large enough in the case of 
polar phase. 



4-1.1. Ferromagnetic phases There are two possible ferromagnetic phases, labelled 
by F and F'. As in F — 1 case, in both of these phases we can use SU(2) instead of 
K x SU{2). 

In .F-phase H F = H F /H F = {I, (-ia z ), (-ia z ) 2 , {-w z f} and thus ni(G/H F ) 
Z4, iT2(G/H F ) = 0. Here a z is the z-component of Pauli matrices and —i<J z describes 
rotation about the z-axis through 180°. Non-trivial vortices are those in which the 
reference spinor rotates through 180°, 360° or 540° about z-axis when the defect line is 
circulated. In SO(3) the isotropy group is {±1}, so G/H = SO(3)/Z2- For a pictorial 
representation of paths in SO(3)/Z2 see [13] . 

In F' phase the order-parameter space is SO(3), and defects are similar to those 
in the ferromagnetic phase of spin-1 condensate. 



4.1.2. Cyclic phase In C-phase we meet an example of a non-commuting first 
homotopy group. A rotation and a gauge transformation of the reference spinor give 

/ e -2-(cos 4 fe-^ + ^sin 2 /3-sin 4 fe^) \ 
e~ ta sin/3( i 



1 



2 £ -i2 7 



£ sin 2 /3-sin 4 fe^) 
- ^ sin 2/3 + sin 2 fe' , 



2 

sin 2 (3 sin 2 7 + ^ (1 + 3 cos 2/3) 



e ia sin /3( sin 2 fe^ 27 + ^ sin 2/3 
e 2ia ( sin 4 f e- 427 + ^ s i n 2 (3 - 



cos 



cos 2 §e i27 ) 

4 §_ e i2j 



) I 



Equating this with | (l, 0, \/2, 0, — l) yields the elements of the isotropy group. H c 
turns out to be a discrete, non-commuting group. Explicitly H c is the union of the 
conjugacy classes shown below. The isotropy group (in U{\) x SO(3)) is isomorphic to 
the tetrahedral group, which is the symmetry group of a tetrahedron. H c is a discrete 
group and thus iri(G/H c ) = H c . Because H c is a non-commuting group we have to 
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use the conjugacy classes of ni(G/H c ) to classify the topologically inequivalent defects 
(5- Two line defects are topologically equivalent if and only if they are characterized 
by the same conjugacy class of the first homotopy group. Defects can still be labelled 
by the elements of the first homotopy group, but if these elements belong to the 
same conjugacy class, corresponding defects can be continuously transformed to one 
another. However, if they belong to different conjugacy classes this is not possible. 
The conjugacy classes are 

<7„(n) = {(n,I)}, Co(n) = {(»,-!)} 

6*2(71) = {(ra, a), (n, -a), (n, b), (n, -b), (n, c), (ra, -c)} 

Cail + n) = {(l + n, d), (§ + n, e), (± + n, /), (| + n,g)} 

C 3 a +n) = {(J + 71, -rf), (| + n, -e), (I + n, -/), (§ + n, -<?)} ^ 

Cf(| + n) = {(I + 7i, d 2 ), (| + 7i, e 2 ), (| + n, / 2 ), (| + 7i, 3 2 )} 

Cf(| + n) = {(I + 71, -d 2 ), (| + 71, -e 2 ), (| + 71, -/ 2 ), (| + 71, - 3 2 )} 

Here n £ Z, a = W(vr,0,0),6 = W(0,7r,f),c = W(0,7r,^),d = W(|,f,f),e = 
W(f,§,±^),/ = W(i^,f,f), 5 = W(f and a 2 = 6 2 = c 2 = d 3 = e 3 = 

/ 3 = ,g 3 = —I. We have also divided the real number part of each group element 
by 2ir. The class Co(n) describes defects in which the phase of the spinor is changed 
by 2im as the defect line is encircled. Notice that only Cq(0) corresponds to trivial 
defects. In the case of Co(n) phase change of 27m is accompanied by a 360° rotation 
about z-axis. For the rest of the conjugacy classes an explicit description of the defects 
is more complicated. For example the element (n, a) in the class 62(71) depicts a defect 
in which the spinor rotates through 180° about z-axis and changes phase by 2mi as 
the line is encircled. Similarly (71, b) € 62(71) describes rotations first through 90° 
about the z-axis and then through 180° about the y-axis together with 2im phase 
change. However, because these defects belong to the same conjugacy class they can 
be continuously transformed into one other. 

The multiplication table of conjugacy classes is shown in table Q] It shows that, 
for example, when we combine defect C2(n) with C 2 (— n) they can either annihilate 
each other (Co(0)) or form defect Cq(0) or 6*2(0), the result depending on how they 
are brought together. 

Defects can be classified further using homology groups ^2lEj- In the presence 
of other line singularities it may be possible to transform two line defects described 
by different conjugacy classes into one another. This is achieved by splitting a defect 
into two parts and combining these beyond a suitable line defect. Elements of ni(M) 
can be grouped into sets in such a way that defects described by elements in the 
same set can be deformed to one another either continuously or by the previously 
described way. The collection of these sets forms a factor group ir 1 (M)/D, where D is 
an invariant subgroup of m(M) generated by elements StS^t^ 1 with 5, r £ wi(M). 
The elements of Ki(M)/D are unions of conjugacy classes. In our case D is the union 
of the conjugacy classes with winding number zero, D = Co(0) U Co(0) U 6*2(0), and 

Tn{G/H c )/D = {Co U Co U C 2 , C7 3 U C^, C 3 2 U C$ }. (5) 

Here we have omitted winding numbers, which are 71, 1/3 + 71 and 2/3 + ?i respectively. 
We see that line defects with the same winding number can be deformed to one another 
either continuously or using a splitting and recombination process. 

From the work of Poenaru and Toulouse ^3] we know that when two line defects 
(described by S, r £ m(M)) cross each other they produce a new line defect connecting 
them. This defect is of the type StS^ 1 ^ 1 . Clearly, if StS^t 1 — 1, line defects can 
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Table 4. The multiplication table of the conjugacy classes of C-phase. Because 
the class multiplication is commutative only half of that is shown. Winding 
numbers have been omitted for clarity. When two classes are multiplied the 
winding number of the resulting class is the sum of the individal winding numbers. 

3 c| cf 
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c 2 


Co 


c 




c 2 


c 2 


6C + 6Cq" + 4C 2 


c 3 


c 3 


3(C 3 + Cj) 




c 3 


3(C 3 + tTJ) 


r 2 

°3 


o 2 

L 3 


3(C| + Cj) 


r 2 

°3 


r 2 

°3 


3(C| + C|) 



3C| + C| 

C| + 3C§ 3C| + C| 

4Uo"+2C 2 4C + 2 C 2 3C3+U3- 

4C + 2C 2 4Co" + 2 C 2 C 3 + 3C3" 3C 3 + Cj" 



Table 5. General forms of the ground state spinors of F = 2 condensate in the 
precence of external magnetic field. Energy is degenerate with respect to angles 
9 and <j>. 

e £ 



Fi e ie (1,0,0, 0,0) 4c-2p 

F 2 e ie (0,0,0, 0, 1) 4c+2p 

F[ e i9 (0,1,0, 0,0) c-p 

F' 2 e ie (0,0,0, 1,0) c+p 

C \ {e W (1 + t), 0, yj - g , 0, e-* C-^) (-! + £)) 

p ^ ( eie / r +^. o . °> 0. ei V 1 - i^a ) ^-3^ 

^ ^ (0. eie V /TT ^' °' ei V 1 " ai^ay • °) 

P e l9 (0, 0,1,0,0) d 



pass through each other without the creation of a new singular defect. In our case 
defects that can be created by making two line defects cross are the trivial defect Co(0) 
and two non-trivial defects, namely Co(0) and 62(0). 

4-2. Non-zero external field 

Ground states were calculated in |H] and are shown in table El However, now 
it should be noted that in the cyclic phase the order-parameter space has a quite 
complicated structure |Hj. Group U(l) x 5*0(2) can act transitively on this order 
parameter space only if the external field is strong enough, and even then there may 
be states which are degenerate in energy but which can not be obtained from the 
reference order-parameter shown in table \^[8\. 

In the ferromagnetic phases and in the Po phase G/H = U{1) and the first and 
second homotopy groups are Z and 0. In C, P and P\ phases G/H — U(l) x U(l) 
and the homotopy groups are Z x Z and 0. Physically this means that we can have 
a vortex in each component of a spinor but only two of them can have independent 
winding numbers. In C-phase, if there are vortices with winding numbers m and n 
say, in the first and third component of the spinor, then there must be a vortex also in 
the fifth component of the spinor. However, its winding number is not free but equal 
to 2n — m. 
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5. Discussions 

In this paper, we have calculated the first and second homotopy groups of the order- 
parameter spaces of spinor condensates with F = 1 and F = 2. The elements of 
these groups correspond to topologically stable singular line and point defects. The 
order-parameter space is identified with the set of degenerate ground state spinors, 
and both non-zero and zero external magnetic field cases are discussed. 

In F = 1 condensate there are two possible phases, ferromagnetic and 
antiferromagnetic. If external field is zero in the former there can be one topologically 
non-trivial line defect but no topologically non-trivial point defects. In the latter 
infinitely many line and point defects, labelled by integers, are possible. 

In F = 2 condensate three different phases, ferromagnetic, polar and cyclic are 
possible. The ferromagnetic phase can be further divided into two phases labelled by 
||(F)||| = 1 or 2. In zero field the former has similar defects to the ferromagnetic phase 
of F = 1 condensate and in the latter there can be three topologically non-trivial line 
defects but point defects are not stable. 

In the absence of an external field the order-parameter space of the cyclic phase 
has a non-commuting first homotopy group. Topologically stable defects are classified 
by the conjugacy classes of this group and are those in which the spinor is suitably 
rotated and its phase changed by an integer multiple of tt/3 as the defect line is 
encircled. Stable point defects are not possible. If external magnetic field is applied 
the symmetry is reduced and non-commutativity of the first homotopy group is lost. 
It also turns out that in the zero field U(l) x 5*0(3) does not act transitively on 
the order-parameter space of the polar phase and thus the defect structure remains 
unsolved. 

For F = 1 and F — 2 condensates, if the external field is non-zero and there is 
only one non-zero component in the spinor, a vortex with an arbitrary integer winding 
number is possible. If there are two or three non-zero components then a vortex in each 
component of the spinor is possible, but only two of these can have an independent 
winding number. In the presence of a magnetic field stable point defects cannot exist. 
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Appendix 

The 3x3 representation matrix corresponding to U(a,P,j) is 



The five dimensional representation matrix is given by (a, 0, 7) = 

cxp(— iaF z ) exp(—i/3F y ) exp(—i~/F z ), where exp(— iaF z ) = diag(e~ l2a , e~ ta , 1, e M , e l2a ) 




(A.l) 
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and 



exp(-i/3F y ) 



& sin 2 /3 





| (cos /3 + cos 2/3) 



|(cos/3 - cos 2/3) 



— sin j3 cos 2 ^ 



sin /3 sin 2 £ 



Y sin 2/3 



^sin 2 /3 
-^psin 2/3 
i(l + 3cos 2/3) 
^ sin 2/3 
^ sin 2 /3 



5 (cos /3 + cos 2/3) 



| (cos /3 - cos 2/3) 



— sin /3 sin 2 ^ 




/ 



(A 
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